Theory of ARPES intensities from the CuC>2 plane 
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We present a theory for the photon energy and polarization dependence of ARPES intensities 
from the CuCh plane in the framework of strong correlation models. We show that for electric field 
vector in the Cu02 plane the 'radiation characteristics' of the O 2p a and Cu 3d x 2_ y 2 orbitals are 
strongly peaked along the Cu02 plane, i.e. most photoelectrons are emitted at grazing angles. This 
suggests that surface states play an important role in the observed ARPES spectra, consistent with 
recent data from Sr2CuCl2 02- We show that a combination of surface state dispersion and Fano 
resonance between surface state and the continuum of LEED-states may produce a precipitous drop 
in the observed photoelectron current as a function of in-plane momentum, which may well mimic 
a Fermi-surface crossing. This effect may explain the simultaneous 'observation' of a hole-like and 
an electron-like Fermi surfaces in Bi2212 at different photon energies. We show that by suitable 
choice of photon polarization one can on one hand 'focus' the radiation characteristics of the in- plane 
orbitals towards the detector and on the other hand make the interference between partial waves 
from different orbitals 'more constructive'. 

71.30.+h,71.10.Fd,71.10.Hf 



I. INTRODUCTION 

Their quasi-2D nature makes cuprate superconductors 
ideal materials for angle resolved photoemission spec- 
troscopy (ARPES) studies, and by now a wealth of ex- 
perimental data is available M. On the other hand, it 
does not seem as if these data are really well-understood, 
the major reason being that we still lack even a rudimen- 
tary understanding of the matrix element effects present 
in these materials. It has recently turned out that ma- 
trix element effects are (or rather: should be) the central 
issue in the discussion of ARPES data. 
ARPES is generally believed to measure the single parti- 
cle spectral function, which near the chemical potential 
fi (and neglecting the finite lifetime) can be written as 



A(fe, w ) = |(* Q p(fe)|c fe)CT |*o)| 2 
0(E QP (k)-ri 6{u 



(E QP (k)-»)) 



Here Eqp denotes the dispersion of the 'quasiparticle 
band', |\&o) an d \*&Qp(k)) denote the ground state and 
quasiparticle state, respectively. In other words, the 
experiment gives a 'peak' whose dispersion follows the 
quasiparticle band Eqp(Jc), with the total intensity of 
the peak being given by the so-called quasiparticle weight 
Z(k) = |(\E'Qp(fc)|c fc , (T |*o)| 2 - For free particles we have 
Z(k) — 1, whence the only reason for a sudden vanish- 
ing of the peak with changing k can be the Fermi factor 
Q(EQp(k) — fi), i.e. the crossing of the quasiparticle band 
through the Fermi energy. Under these circumstances, 
it would be very easy to infer the Fermi surface geome- 
try from the measured photoelectron spectra, and indeed 
this very assumption, namely that a sudden drop of the 
photoemission intensity automatically indicates a Fermi 
level crossing, has long been made in the interpretation 



of all experimental spectra on metallic cuprates. 
Several experimental findings have shown, however, that 
this assumption is not tenable in the cuprates. The first 
indication comes from the study of the insulating com- 
pounds Sr 2 CuCl 2 2 § and Ca 2 Cu0 2 Cl 2 §]. Although 
these insulators cannot have any Fermi surface in the 
usual sense, which means that the factor 0(EQp(k) — fi) 
is always equal to unity, the experiments show that also 
in these compounds the quasiparticle peak disappears as 
one passes from inside the noninteracting Fermi surface 
to outside. Thereby a particularly striking feature of the 
experimental data is the sharpness of the drop in spectral 
weight ||, which is for example along (1, 1), quite com- 
parable to the drops seen at the 'Fermi level crossings' in 
the metallic compounds. The only possible explanation 
for this phenomenon is a quite dramatic fc-dependence of 
the quasiparticle weight, Z(k). Apparently in these com- 
pounds we have very nearly Z(k) oc ®{Ef ree {k) — fi), i.e. 
the k dependence of Z(k) resembles that of the nonin- 
teracting system. 

This result immediately raises the question as to how sig- 
nificant the 'Fermi level crossings' observed in the metal- 
lic compounds really are. That they may, in some cases, 
have little or no significance for the true Fermi surface 
topology has been demonstrated by the recent contro- 
versy as to whether the Fermi surface in the most exhaus- 
tively studied compound, Bi 2 Sr 2 CaCu 2 0s+,5, is hole-like 
(as inferred from a large number of studies |(| with 
photon energy 22eV) or electron-like (as concluded by 
several recent studies ]7|-pT| at photon-energy 33eV). 
Here it should be noted that the true Fermi surface topol- 
ogy is an intrinsic property of the material which can 
under no circumstances change with the photon energy. 
It follows from these considerations that to extract any 
meaningful information from angle-resolved photoemis- 
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sion we need an understanding of the matrix elements 
and their fc-dependence, as well as other effects which 
might possibly influence the intensity of the ARPES sig- 
nal. Motivated by these considerations, we have per- 
formed a theoretical analysis of the spectral weight of 
strongly correlated electron models. In section II we de- 
rive a simple expression for the photoelectron current, 
which can be applied e.g. in numerical calculations for 
strong correlation models. In section III we specialize 
this to the CuC>2 plane, in section IV we discuss the an- 
gular 'radiation characteristics' of a Zhang-Rice singlet 
and how these could be exploited to optimize the pho- 
toemission intensity from the respective states. We also 
show that photoelectrons are emitted predominantly at 
small angles with respect to the CuC-2 plane. In sec- 
tion V we point out that this may lead to the injection of 
these photoelectrons into surface resonance states. Fano- 
resonance between the surface resonance and the contin- 
uum of LEED states then leads to a strong energy de- 
pendence of the ARPES signal and we show that already 
a very simple free-electron model can explain the exper- 
imental energy dependence of the first ionization states 
in Sr 2 CuCl202 measured recently by by Diirr et al. |l|] 
surprisingly well. In section VI we describe how the inter- 
play between surface state dispersion and Fano-resonance 
between the processes of direct emission and emission via 
a surface state can mimic a Fermi level crossing where 
none exists, and suggest that the apparent change in 
Fermi surface topology with photon energy may be due 
to such 'apparent Fermi surfaces'. Section VII contains 
our conclusions. 

II. PHOTOEMISSION INTENSITIES FOR 
STRONG CORRELATION MODELS 

Photoelectrons with a kinetic energy in the range 
10— lOOeV have wavelengths comparable to the distances 
between individual atoms in the CuC>2 plane. It follows 
that whereas the eigenvalue spectrum of the plane prob- 
ably can be described well by an effective single band 
model |p^| , this is not possible for the matrix elements. 
We necessarily have to discuss (at least) the full three- 
band model. 

Our first goal therefore is to derive a representation of 
the photoemission process in terms of the electron anni- 
hilation operators for the Cu 3d x 2_ y 2 and O 2pa orbitals 
in the CuC>2 plane. In other words, we seek an operator 
of the form 

Ck.a = m(d x i^ y i)d k ^ + m(p x )p x ,h,<7 + m(p y )p v , h ,a 

such that the single particle spectral density of this op- 
erator 



evaluated for the correlated electron model in question 
reproduces the experimental photoemission intensities. 
Here p x .k.a an d P y ,k.a an d arc the annihilation operators 
for an electron in an x- and y-dirccted cr-bonding oxygen 
orbital, dk,a annihilated an electron in the d x 2_ y 2 orbital. 
To that end, let us first consider the problem of a single 
atom (which may be either Cu or O). The calculation is 
similar as outlined in Refs. jl3| and iQ. We want to 
study photoionization, i.e. an optical transition from a 
localized valence orbital into a scattering state with en- 
ergy E. The dipole matrix element for light polarized 
along the unit vector e reads: 

m, w = J dr $* f (r) (e • r) %(r) 

Here the initial state is taken to be a CEF state with 
angular momentum V = 1,2 and crystal- field label a — 

PxiPy ) d x 2_y2 .... 

**( r ) = ^ J C<xm'Yv t m'{r ) Rn V {r). 
m' 

We consider this state to be 'localized', that means the 
radial wave function R n i (r) is zero outside the atomic ra- 
dius tq. For the wave function of the final (^scattering) 
state we choose: 

{ Y lm (r°) 2sin(fcr + <y,-^) r > r 
*/(r) = { 

{ Y lm (r°) ±ReM r < r 

The real functions R n i>(r) and Ri(E) both are a solution 
of the radial Schrodinger equation for a suitably chosen 
atomic potential. The scattering phase Si and the pref- 
actor — are determined by the condition that the wave 
function $ At) and its derivative be continuous at r = r Q . 
Details are given in Appendix I. We note that the scat- 
tering phase Si also plays an important role in the inter- 
pretation of EXAFS spectra (where it is usually called 
the central atom phase shift) and thus could in principle 
be determined experimentally (although at the relatively 
high photoelectron kinetic energies important for EXAFS 
are not ideal for ARPES). 
Representing the dipole operator as 

we can rewrite the dipole matrix element as 

i 

dl, m = Cam ' C 1 ^, m; 2', m') y*, m -m'( e )- (!) 

m' 

Here the radial integral is given by 



2 



i r° 

Ri,i>{E) = - / drr 3 Re,i{t) R n v{r) 
n Jo 

and the following abbreviation for the angular integrals 
of three spherical harmonics has been introduced: 

J dnYC m (r )Y 1 ^r a )Y l , m ,(r°) = 
j 3 

S m , m '+n \ — c (l,m\ l',m'). 

V 47T 

The constants c 1 (Z, m; V , m!) are well-known in the the- 
ory of atomic multiplets and tabulated for example in 
Slater's book p5[ . Knowing the atomic potential the ra- 
dial integral and thus the entire matrix element could 
now in principle be calculated. 

Next, we assume that the radial part of ^ /(r) at large 
distances is decomposed into outgoing and incoming 
spherical waves. Observing the outgoing spherical wave 
at large distance under a direction defined by the polar 
angles @k and it may locally be approximated by a 
phase shifted plane wave: 



¥/(r) 



2mE o 



where k - 

Let us next consider an array of identical atoms in the 
(x, y)-plane of our coordinate system (which we take to 
coincide with the CuC>2 plane in all that follows) . To de- 
scribe the final state after ejection of an electron from the 
atom at site j we would simply have to replace through- 
out r — > r — Rj in the above calculation, where Rj 
denotes the position of the atom. If we want to give 
the created photoholc a definite in-plane momentum fen, 
however, we have to form a coherent superposition of 
such states, that means weighted by the Bloch factors 
_i e -ik r Rj w h ere denotes the number of atoms in 

VN 

the plane. At the remote distance r we consequently re- 
place r — > r — Rj, leaving the polar angles 8^ and 
unchanged, multiply by ~ afc n Rj and sum over j. 
The photoelectron wave function at r then becomes 



(k — fey )-Rj 



F/,™(fe u ) H) (+i e 



NSk tl +G„,k x 



(2) 



Here k x y denotes the projection of k onto the (x, y)- 
plane and Gn is a 2D reciprocal lattice vector. An im- 
portant feature of this result is the fact that by creating a 
photohole with momentum kn (which must belong to the 



first Brillouin zone) the photoelectrons may well be emit- 
ted with the 3D momentum (fey + G\\,k±), that means 
the parallel momentum component of the photoelectrons 
need not be equal to fey. 

Summing over all possible partial waves (l,m) and in- 
troducing the abbreviation Riy(E) = Ri^(E)e lSl the 
electron current per solid angle at r due to the in-plane 
orbitals of the type (I', a) finally becomes 



ANhk 



J = 



d lm R u ,(E)Y lm (k a ) H) i+1 | 2 



Here N should be taken equal to the number of unit cells 
within the sample area illuminated by the photon beam. 
After some algebra (Appendix II) this can be brought to 
the form 



J = — (vi, a (k°)-e) 
i 



in 



(3) 



Here all the angular dependence on the shape of the orig- 
inal CEF-level has been collected in the vectors Vi :Ct - it 
is important to note, that these vectors are obtained by 
standard angular-momentum recoupling so that the re- 
sulting angular dependence is exact. It is only the 'radial 
matrix elements' Riji (E) which have to be calculated ap- 
proximately and thus may be prone to inaccuracies. 
So far we have limited ourselves to one specific type of 
orbital in the plane. If we allow photoemission from any 
type of orbital we simply have to add up the prefactors 
of the plane-wave states before squaring to compute the 
photocurrent. 

Summarizing the preceding discussion we can conclude 
that the proper electron annihilation operator to describe 
the photoemission process would be 

c fe|| , CT = J2e- ik " Ra (Ri lV(a) {E) K Q (fc°) • e) ) c a , kh<7 

l,a 

= y^v(o!) ■ e c a ^ u <j- (4) 



Here R a denotes the position of the orbital a within 
the unit cell. Suppressing the spin and momentum in- 
dex for simplicity, the total ARPES intensity then would 
be given (up to an energy and momentum independent 
prefactor) by 

I(u) = (v(«) ' e)* (v(/3) • e) QR a ,^ - *«+) (5) 

a, 13 



RaA z ) = -<o|4 



w - (H - Eq) - i0+ 



c p \0). 



(6) 



Whereas the scalar products (v(a) • e) take into account 
the interplay between the real-space shape of the orbitals, 
the polarization of the incident light and the direction of 
electron emission, the spectral densities R a .p(z) incorpo- 
rate the possible many-body effects in the Cu02 planes - 
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we thus have the desired recipe for studying photoemis- 
sion intensities in the framework of a strong correlation 
model. 

To conclude this section we note that we have actually 
performed only the first stage of the calculation within 
the so-called three-step model. We give a brief list of 
complications that we have neglected: the emission from 
lower planes than the first and the extinction of the re- 
spective photoelectron intensity, any diffraction of the 
outgoing electron wave function from the surrounding 
atoms, the refraction of the photoelectrons as they pass 
the potential step at the surface of the solid. It is thus 
quite obvious that our theory is strongly simplified. 

III. APPLICATION TO THE THREE-BAND 
MODEL 

We now want to discuss some consequences of the 
results in the preceding section, thereby using mainly 

( °' 

H= °' 

-2it pd sm(^f), 

\ -2ii sp sin(Y) 



the standard three-band Hubbard-model to describe the 
CuC>2 plane. We first consider the noninteracting limit 
U = 0. Here we use a 4-band model which was intro- 
duced by Andersen et al. jl6) to describe the LDA band- 
structure of YBa2Cu304. In addition to the Cu id x 2_ y 2 
orbital and the to er-bonding O 2p orbitals this model in- 
cludes a Cu 4s orbital, which produces the so-called £' and 
t" terms in the single-band model - which are essential to 
obtain the correct Fermi surface topology. Suppressing 
the spin index the creation operators for single electron 
eigenstates read 

lk,u = a h»Pk,x + a %»A, y + "3,^4 + a i<u s{ 
where a.^ denotes the v th eigenvector of the matrix 



0, 2it pd sin(^) 2^^ s^) sin(^)\ 

0, -2it pd sm(^) 2it sp sin(^) 

2it pd sm(^) A 

-2it sp sin(^f) A s J 



surface normal (as will be shown below). 



15 30 45 60 75 90 



Values of the parameters t pd t sp and A are given in d x 2_ y 2 orbital produce virtually no intensity close to the 
Rcf. jjj). 

If we consider the planar momentum fcy = G + fe, where 
k is in the first BZ and G is a reciprocal lattice vector, 
the photoemission intensity for the v th band becomes 

I oc \[v[p x ) ■ e) e- lG */ 2 al + (v(p y ) ■ e) e^ G ^ 2 a* 2 
+ {v(d x 2_ y 2) ■ e) a* 3 + (v(s) ■ e) a* 4 \ 2 . (7) 

Using this we proceed to a detailed comparison with 
the work of Bansil and Lindroos p7| . These authors 
have performed an extensive first-principles study of the 
ARPES intensities in Bi2212, thereby using the more re- 
alistic one-step model of photoemission and a complete 
surface band-structure, both for the initial and final state 
wave function. Amongst others, Bansil and Lindroos 
studied the variation of the peak- weight along the Fermi 
surface for given direction of the polarization vector e. 
Their results are shown in Figure [j and compared to the 
above theory. Obviously the overall trends seen by Ban- 
sil and Lindroos are reproduced reasonably well by the 
present theory and we want to give a brief discussion of 
the mechanisms which lead to these trends. 
To begin with, it is to simplest approximation the 'oxy- 
gen content' of the wave function, which determines the 
ARPES intensity. This is not so much due to the small- 
ness of the radial matrix elements for Cu (which are quite 
comparable to the ones for oxygen, see the Appendix), 
but rather the fact that the partial waves emitted by a 





FIG. 1. ARPES intensity at the Fermi energy (in ar- 
bitrary units) versus Fermi surface angle $ (see lower right 
figure), for three different polarization directions. The curves 
were computed from the tight-binding model (solid line) and 
ZRS (dashed line). The photon energy is hv = 22 eV and the 
curves are normalized to $ = in panel the upper left panel. 
The points represent the integrated intensities from first the 
first principle calculation of Bansil and Lindroos The 
lower right panel shows the geometrical details. 
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For an electron energy of « 20 eV the detector would 
have to be placed at w 20° degrees from the surface nor- 
mal (neglecting the refraction by the potential step at 
the surface), and the d x 2_ y 2 orbital emits practically no 
electrons into this direction. 

Next we note that light polarized along the (1, 0)direction 
can only excite electrons from p^-type orbitals, light po- 
larized along the (1,1) direction will excite the combina- 
tion p x +p y , whereas light polarized along (1,-1) excites 
Px — Py Finally one has to bear in mind that near (tt, 0) 
the tight-binding wave function contains only p^-type or- 
bitals (the mixing between p v (tt, 0) and d x 2_ y 2(ir,Q) be- 
ing exactly zero). Therefore the polarizations (1,1) and 
(1,-1) have equal intensity, whereas (1,0) gives maxi- 
mum intensity. Near (0, tt) in the other hand, the tight- 
binding wave function contains only p y , whence the po- 
larizations (1,1) and (1,-1) again have equal intensity, 
whereas (1,0) this time gives no intensity. For k || (1,1) 
exciting p x + p y gives no intensity, because this combi- 
nation does not mix with d x 2_ y 2 and hence is not con- 
tained in the tight-binding wave function, whereas excit- 
ing p x — p y gives high intensity. The intensity for po- 
larization (1,0) at these momenta is approximately 1/2 
of that for (ir,0). These simple considerations obviously 
explain the overall shape of the curves in Figure [l] quite 
well. The main discrepancy between the Bansil-Lindroos 
theory and our calculation is the behaviour of the curve 
for polarization (1,-1). Actually, this is the only curve 
which is not determined by symmetry alone, and sub- 
tle details of the wave function become important. We 
believe that the parameterization of Andersen et al. |T(| 
does give the correct dispersion, but not necessarily the 
right wave functions. Next we consider the intensity vari- 
ation along various high-symmetry lines in the Brillouin 
zone, shown in Figure |2|. 



Bansil and Lindroos, 22eV a 

ZRS, 22eV (8eV LP.) 

TB, 22eV (8eV LP.) 

CPT, 22eV (8 eV LP.) 




(0,0) (n,0) (Jt,lt) (0,0) 



k 

FIG. 2. Spectral weight along (0,0) -> (tt, 0) and 
(0, 0) — > (tt, tt) compared to first principles calculation fH] ] 
(triangles). ZRS (solid line), tight-binding (long dashes) and 
CPT (short dashes) method show qualitatively good agree- 
ment. The strong decrease of the spectral weight in CPT 
towards k = (tt, 0) may result from high doping (S = 25%). 




FIG. 3. Orbitals used for constructing a ZRS and signs of 
the different lobes. 

Again, there is reasonable agreement between Bansil and 
Lindroos and our theory. All in all the comparison shows 
that also in the results of Bansil and Lindroos the be- 
haviour of the intensity is to a considerable extent de- 
termined by the 'radiation characteristics' of the orbitals 
in the Cu02 plane and the relative weight of the p x and 
p y orbitals. Additional complications like multiple scat- 
tering of the photoelectrons, Bragg scattering from the 
BiO surface layer etc. do not seem to have a very strong 
impact on the intensity variations, not even at the rela- 
tively low photon energy of 22 eV. Despite its simplic- 
ity we therefore believe that our theory has some merit, 
particularly so because it allows (unlike the single par- 
ticle calculation of Bansil and Lindroos) to incorporate 
the effects of strong correlations. One weak point of all 
single-particle-like calculations for the Cu02-plane is the 
following: since (in electron language) the band which 
forms the Fermi surface is the topmost one obtained by 
mixing the energetically higher Cu 3d x 2_ y 2-orbital with 
the energetically lower O 2pcr orbitals it is clear, that the 
respective wave functions have predominant Cu 3d x 2_ y 2 
character. This is exactly opposite to the actual situa- 
tion in the cuprates, where the fist ionization states in the 
doped and undoped case are known to have predominant 
O 2p character. We now consider the case of large U. 
An exact calculation of the single-particle spectra (||) is 
no longer possible in this case and we have to use various 
approximations. First we study an isolated Zhang-Rice 
singlet (ZRS) in a single Cu04 plaquette, see Figure ||. 
The bonding combination of O 2pcr-orbitals is 

P{a = -P3,a+Pl,a) 

and the single-hole basis states (relevant at half-filling) 
can be written as 

|1) = Pb : a\f ull )> 

\2) = d„\full). (8) 

Here \full) denotes the Cu3d 10 <g) 4 02p 6 state. The sin- 
gle hole ground state then is ^) = a\l) + /3\2) where 
(a, [3) is the normalized ground state eigenvector of the 
matrix 



5 



Two-hole states are obtained by starting from the basis 
states 

I 1 ') = Pb,tPb,l\full) 

\2') = ^=(p b , 1 d L + dfpb^lfull) 

\3') = d^full) (10) 

The two-hole ground state of the plaquette reads 
l^)) = a '\V) + (3'\2')+i\3'), where (a',/?', 7') is an 
eigenvector of 

/ 2V2t pd \ 
H = 2V2t pd -A 2V2t pd . (11) 
V 2V2t pd -2A + U J 

If we want to make contact with the t — J model, where a 
'hole' at site i stands for a ZRS in the plaquette centered 
on the copper site i, we have to incorporate a phase fac- 
tor of e _ * fe H' Ri , into the definition of | \E r Q 2/l ' ) ) , where Ri 
is the position of the central Cu orbital. In other words, 
the matrix element for the creation of a 'hole' in the t — J 
model is 

m ZRS = e^li'^^^lc^l^). (12) 

The matrix elements for the creation of e~ ^' R ^ p b>cr \ full) 
and e~ lh \y R id a \full) are 

rn b = i (v(p y ) sin(y ) - v( Px ) sin( y ) j • e, 

m d = v(d x 2_ y 2) ■ e. (13) 

Finally, the matrix clement for creation of a ZRS from 
the single-hole ground state becomes 

f3'* 

mzRS = a'*am b + —=(am d + f3m b ) 
\2 

+l'*Pm d . 

Using this expression, the intensity expected for a ZRS 
with momentum k can be calculated as a function of 
photon polarization and energy. This will be discussed 
in the next sections. 
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Ocu 3d x 2. 
O O 2p Xjy 



FIG. 4. Three band Hubbard cluster used in the Lanczos 
diagonalization. The hopping across the cluster boundary 
(dashed line) is treated to lowest order strong coupling per- 
turbation theory. 



In comparing the intensities calculated from our above 
theory for a single ZRS to experiment one would be im- 
plicitly assuming that the first ionization states can be 
described completely by a coherent superposition of ZRS 
in a single plaquette. This need not be the case and in 
order to get at least a rough feeling for the effects of 
'embedding the plaquette in a lattice' we use the tech- 
nique of cluster perturbation theory (CPT) to study the 
extended system. CPT, which was first suggested by 
Senechal et al. is a technique to 'extrapolate' the 
single-electron's Green's function calculated on a finite 
cluster to an infinite periodic system. It is based on an 
perturbative treatment of the intercluster hopping. The 
cluster we used for the present calculation is a quadratic 
arrangement of four unit cells each containing one Cu 
d x 2_ y 2 and two O p x ^ y orbitals as shown in figure f§ The 
Green's function 



($o|c i(T - 



1 



cL|$ > + (*o|c 



1 



ia z-H 



Cia |$ ) 



of this cluster with open boundary conditions is calcu- 
lated by the Lanczos method. The boundary orbitals 
of the cluster are connected to the adjacent cluster by 
the Fourier transform of the intercluster hopping V(Q), 
where Q is a "superlattice" wave- vector restricted to the 
smaller Brillouin zone formed by the now enlarged lat- 
tice of the clusters. The perturbative treatment of the 
intercluster hopping in lowest order yields an RPA like 
expression for the approximate Green's function. 



Qija (Q, 



1 



G~i(uj) a -V(Q) 



This mixed representation can then be transformed into 
momentum space by a Fourier transformation thereby 
taking into account the geometry of the cluster: 



Gcpt 1 



This result is exact for U = and has been shown [Q to 
produce good results also for large and intermediate U . 
In the following sections we will use this technique for 
the approximate calculation of the single-particle spec- 
tral densities as a valuable cross-check for the calcula- 
tions based on a single ZRS. While CPT still is far from 
being rigorous, it includes some effects which result from 
embedding the ZRS in a lattice and, as we will in fact 
see, it always gives results which are very similar to those 
for a ZRS. 

To conclude this section we want to give a simple esti- 
mate for the photoelectron kinetic energy T as a function 



G 



of photon energy hv. It is easy to see that in the absence 
of any potential step at the surface of the solid we would 
have 



T = hv + E, 



(N) 




E 



(N-l) 



(14) 



where Eq N ^ is the energy of the ground state of the solid 



with N electrons and E, 
We approximate 



(N-i) 



the final state of the solid. 



E 



(N) 



E, 



(N-l) 
k 



(Tplh Tp2h\ w 
1-^0 _ ^0 ) ~ h 2 P 



2h \ 



(15) 



Here Eft h are the ground state energies of a CuC>4 pla- 
quette with n holes, obtained by diagonalizing the ma- 
trices (||) and (11). In these matrices the energy of the 
02p level, Ei p , has been chosen as the zero of energy, 
whence it has to be taken into account separately. For 
the 'standard values' of the parameters t p d — — 1.3 eV, 
A = 3.6 eV and U = 10.5 eV this gives Eft 11 - Ef 1 = 
1.93 eV. Using the estimate Ei v w — 10 eV from our 
atomic LDA calculation we find 



T w hv - 8 eV, 
which we will use in all that follows. 



(16) 



Mesot et. al, 22eV 
Mesot et. al., 34eV 
Bansil and Lindroos , 22eV 
CPT, 22eV 
CPT. 34eV 
TB: 22eV 




expressions for the vectors vi a from Appendix II with 
the 'form factor' of a ZRS, (|l3"|), the expression for the 
photocurrent can be brought to the form 



(17) 



i.i' 



We remember that fey is the momentum of the photo- 
hole (which is within the first Brillouin zone of the CuC>2 
plane) , whereas k is the momentum of the escaping pho- 
toelectron. Using the expressions from Appendix II it is 
now a matter of straightforward algebra to derive the fol- 
lowing expressions for k\\ = (k, 0), and the 'optimal' pho- 
ton polarization for momenta along (1,0), e = (1,0,0): 



«1,2 ' e 



i> 2 ,i ' e = 



V3,2 ■ e 



sin(fe/2) 



m/— ain(e), 



sin(fe/2) /3 ,„^, 1 




M/^sin(9) ( ^cos(2e)-i; 



(18) 

(19) 
(20) 

(21) 



Figure [6] shows the 0-dependence of these expressions. 
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FIG. 5. Spectral weight along (0, 0) — » (it, 0) and 
(0, 0) — ► (n, 7r) compared to experimental data from refer- 
ences 

As a first application, Figure ^Jthen shows the intensity of 
the topmost ARPES peak calculated by CPT and com- 
pares this to experimental data. Reasonable qualitative 
agreement can be found, especially for the CPT calcu- 
lation, although the shift of the spectral weight maxi- 
mum towards k — (tt, 0) for 22 eV is not reproduced. 
The much better agreement for E p hoton = 34eU suggests 
that for the lower photon energy of 22 eV additional 
effects (such as multiple scattering corrections or Bragg- 
scattering from the BiO top-layer) are more important. 



IV. APPLICATION TO EXPERIMENT 

Coming back to the theory for the ZRS we can already 
draw conclusions of some importance. By combining the 



FIG. 6. 'Radiation characteristics' \viy ■ e| 2 for the differ- 
ent partial waves from a ZRS with momentum (it,0). The 
figure shows the O-dependence of the respective partial wave 
within in the x — z plane. The polarization vector e of the 
exciting light is assumed to be in a;-direction. 

Obviously a ZRS with momentum near (tt, 0) emits pho- 
toelectrons predominantly parallel to the Cu02 plane if 
it is excited with light polarized along (1,0) within the 
CuC>2 plane. The situation is similar, though not as pro- 
nounced, for other momenta, such as (§,§)• It should 
be noted, that the vectors Vij* (fey , k) are computed ex- 
actly namely by elementary angular-momentum recou- 
pling. The fact that a ZRS with momentum near (tt, 0) 
emits photoelectrons predominantly at small angles with 
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respect to the Cu02 plane therefore is a rigorous result. 
The radial matrix elements R^i' (E) which also enter in 
( |l7| ) actually tend to suppress the emission close to the 
surface normal even more. Namely one finds (Appendix 
I) that i?o,i ~ 0.2 i?2,i, that means the s-like partial 
wave (which would contribute strongly to emission at 
near perpendicular directions, see Figure ^ has a very 
small weight due to the radial matrix elements. We note 
in passing that the smallness of the ratio i?o, 1/^2,1 is 
well-known in the EXAFS literature pE ] . 
Combining the partial waves in Figure! with the proper 
radial matrix elements Ri } i> (E) as in ( p.7[ ) we expect to 
obtain a curve with a minimum for some finite O (mainly 
due to to the node in the dominant d-like / = 2 partial 
wave emitted by the 2p /' = 1 orbital (see Figure ^|) . This 
may in fact explain a well-known fij effect in ARPES, 
namely the the relatively strong asymmetry with respect 
to (it, 0) of the ARPES intensity for kn along (1, 0) which 




(0,0) 



(it,0) 



(2ic,0) 



(3n,0) 



FIG. 7. ARPES intensity as a function of momentum along 
(1,0) for different photon energies, calculated by for an iso- 
lated ZRS. The inset shows the radiation characteristics of a 
ZRS at (71-, 0) and the directions at which a detector would 
have to be placed at photon energy 22 eV to observe the re- 
spective momentum. Thereby the work function is neglected. 

is seen at 22 eV photon energy but not at 34 eV. Figure 
shows the calculated intensities as a function of in-plane 
momentum. For hv k, 22 eV the anglesG which would 
be appropriate to observe fey slightly beyond (jr, 0) (i.e. 
in the second zone) are such, that one is 'looking into the 
node' of the radiation characteristics of the ZRS, whereas 
for smaller fen one is looking at the maximum (see the in- 
set). Hence there is a strong asymmetry around (it, 0). 
Increasing the energy to > 30eV the angles 9 becomes 
smaller, one is no longer sampling the node and the the 
intensity is much more symmetric around (7r, 0). 
As shown in the preceding discussion, our theory repro- 
duced, despite its simplicity, some experimental features 
seen in the cuprates. We therefore proceed to address 
potential applications in experiment. Thereby the main 
goal is to find experimental conditions under which the 
ZRS-derived state at a given momentum kii can be ob- 
served with the highest intensity. To that end, we can 
vary different experimental parameters, mainly the direc- 
tion of the photon polarization and the Brillouin zone in 



which we are measuring. 

We first consider the case of normal incidence of the light, 
that means the electric field vector e is in the CuC-2 
plane. Rotating e in the plane will change the inten- 
sity from the ZRS-derived states in a systematic way, 
and for some angle <&e it will be maximum. Along high 
symmetry directions like (1, 1) or (1,0) the optimal po- 
larization can be deduced by symmetry considerations, 
because the ZRS has a definite parity under reflections 
by these directions. For momenta k along (1,1) e has 
to be perpendicular to fc, whereas along (1, 0) it must be 
parallel H. For nonsymmetric momenta, however, this 
symmetry analysis is not possible and one has to calcu- 
late the optimum direction. We note that apart from 
merely enhancing the intensity of the ZRS, knowledge of 
the polarization dependence of the ARPES intensity of 
the 'ideal' ZRS would allow also to separate the original 
signal from ZRS-derived states from any 'background', 
that means photoelectrons which have undergone inelas- 
tic scattering on their way to the analyzer. Since it is 
plausible that these background electrons have more or 
less lost the information about the polarization of the 
incoming light, their contribution should be insensitive 
to polarization. Taking the spectrum at the 'optimal 
angle' for the ZRS and perpendicular to it then would 
allow to remove the background by merely subtracting 
the two spectra (provided one can measure the absolute 
intensity). Along the high symmetry directions (1,0) and 
(1,1) the feasibility of this procedure has recently been 
demonstrated by Manzke et al. [^0| and using the cal- 
culated optimal polarizations this analysis could be ex- 
tended to any point in the BZ. 

For the special case of the ZRS, and neglecting the con- 
tribution of copper altogether, it is possible to give a 
rather simple expression for the optimal angle. The ma- 
trix element for creating the bonding combination of O 
2pcr orbitals is 

/ }~ k \ 

mzRS oc i ( sm($ E ) sin(^) - cos($ B ) sin(y ) J . (22) 

We thus find the angles which give minimum and maxi- 
mum intensity: 

k k 
sin($ Bjmin ) sin(-^) - cos($ min ) sin(y) = 

k k 
cos($ E , max ) sin(^) + sin($ max ) sin(y ) = (23) 



We see that always <&E,min = ®E,max + f ■ Moreover, 
in this approximation the intensity is exactly zero for 
&E,mim and the angle for maximum intensity is 



E.max 



= —arctan 



sin(f ) 
sin(4f) 



(24) 



This expression takes a particularly simple form along 
the line (it, 0) — > (0,7r), where $ max = -k y /2. This 
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reproduces the known values <E> max = along (1,0), 
$ma X = ~f along (1,1) and $ max = -f along (0,1). 
Despite its simplicity formula ( p4| ) gives a quite good es- 
timate for the optimum polarization. Numerical evalua- 
tion shows that the lines of constant $ max in fc-space are 
to very good approximation straight lines through the 
center of the Brillouin zone. This means that in practice 
one could adjust the polarization angle $_e once to ci- 
ther $>E,max or <&E,min, and then scan an entire straight 
line through the origin by varying the emission angle Qk 
without having to change the polarization along the way. 
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FIG. 8. Optimal in-plane polarization for observing the 
ZRS at the respective position in the Brillouin zone by 
ARPES with photoelectron kinetic energy energy. The calcu- 
lation is done for an isolated ZRS. 

Figures || and ^| then show the optimal angle for observa- 
tion of the first ionization state for momenta in the first 
Brillouin zone. For perpendicular polarization the inten- 
sity practically vanishes. In the Figure we compare the 
'isolated ZRS' and the 'embedded ZRS' whose spectra 
are obtained by CPT. Both methods of calculation give 
very similar results for the optimal angle, which in turn 
agrees very well with the simple estimate |24|. Since along 
the high symmetry lines (1, 0) and (1, 1) the optimal $e 
is determined by symmetry alone there is obviously not 
so much freedom to 'interpolate' smoothly between these 
values. 

Next, moving to a higher Brillouin zone allows to enhance 
the intensity of the ZRS, as can be seen from Table | 
and y. Table | shows the ratio I/Iq of the intensity ob- 
tainable by measuring the ZRS at (f , f) at the 'original 
position' (there the polarization has to be perpendicu- 
lar to (1,1) by symmetry) and the intensity that can be 
obtained by measuring in a higher zone with optimized 
polarization, both within in the CuC>2 plane and, for later 
reference, also out of plane. 
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FIG. 9. Optimal in-plane polarization for observing the 
first ionization state of the correlated Cu02 plane. The in- 
tensities are calculated by integrating the spectral weight ob- 
tained by CPT within 300meV. 

The same information is given for (tt, 0) in table || Ob- 
viously, an enhancement of a factor of 2 or more can 
be achieved by proper choice of the experimental condi- 
tions. In view of the large 'background' in the experi- 
mental spectra, even such a moderate enhancement may 
be quite important. 

We note that the physical origin of the variation of in- 
tensity with $ e and order of the Brillouin zone 
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FIG. 10. 'Radiation characteristics' within the (x—z) plane 
of a ZRS excited with light polarized in the {x — z) plane. The 
polar angles are Ob = ^ (dashed line) and Ob = Qs.opt ~ f 
(full line) . The skew line shows the direction where the detec- 
tor would have to be placed for the photon energy of 22 eV 
(work function neglected). 
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is the interference between the processes of creating a 
photohole in p^-like and Py-like O 2pa orbitals. By 
changing <£>£ the relative intensity and phase between 
photoholes in the two types of orbitals can be tuned. 
Thus, one can always find an angle $>E,max where the in- 
terference between these two types of orbital is maximally 
constructive. This angle depends on the relative phase 
between the two orbitals in the wave function and thus is 
specific for the state in question. Next, we study a quite 
different effect, namely what happens if we tilt the polar- 



ization vector e out of the CuC>2 plane. To illustrate the 
usefulness of this, we again consider a momentum along 
the high symmetry line (0, 0) — > (tt, 0) and assume that 
the polarization vector e is within in x — z plane It fol- 
lows from symmetry considerations that any component 
of the light perpendicular to this plane ('s-polarization') 
cannot excite any ZRS-derived states. In other words we 
assume that but that @e is variable. The contri- 
bution from the bonding combination of O 2pa orbitals 
then becomes 



m Z RS = -«6 



R.2id 3z 2_ r 2 —R,2id x 2_ l 

15 V5 



sin(6) + —=R 2 id xz cos(0) 
V5 

1 - 



1- 1 ~ / 3 ~ \ 

-—^R32lPx ~ -^=i?323/x(52 2 -l) + y 1^R323 f x^ -3xy 2 J SHl(6) + ~^R323fz(x 2 -y 2 ) C0s(6) 



^ and 6 = i 



P+^f. Tilting 



where £i = a * 
the electric field vector out of the plane thus admixes the 
d xz harmonic into the radiation characteristics, which has 
its maximum intensity at an angle of 45° with respect to 
the CuC>2 plane. Clearly, this enhances the intensity at 
intermediate 0&, particularly so if constructive interfer- 
ence with the d x 2_ y 2 and d 3z 2_ r 2 harmonics (which are 



excited by the x-component of e) occurs 
on one hand and the cL,2_„2 and d 



Since the d xz 
y * uiiu ^ 3z 2_ r 2 on the other 
have opposite parity under reflection by the y — z plane, 
constructive interference for some momentum (fc, 0) auto- 
matically implies destructive interference for momentum 
(— A;,0). Tilting the electric field out of the plane thus 
amounts to 'focusing' the photoelectrons towards the de- 
tector. This effect is illustrated in Figure |l^, which shows 
the angular variation within the (x — z) plane of the pho- 
tocurrent radiated by a ZRS. 
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FIG. 11. Intensity of ZRS derived states with different mo- 
menta along the (1, 0)-direction, as a function of the polar an- 
gle Qe- The light is polarized in the (x — z) plane, i.e. <E>e = 0, 
the photon energy is 22 eV (work function neglected). 



Figure |Tl| then shows the intensity as a function of Qe 
for various momenta along (0,0) — > (tt, 0). Remarkably 
enough, the enhancement of current towards the detector 
may be up to a factor of 3 as compared to polarization in 
the plane. Exploiting this effect could enhance the pho- 
toemission intensity considerably in the region around 
(tt, 0) - which might be important, because this is pre- 
cisely the 'controversial' region in fc-space which makes 
the difference between hole-like [[l]j4]-|(| and electron-like 
[ 7|— 10|| Fermi surface in Bi2212 and where bilayer-splitting 
[21p2[| should be observable. One caveat is the fact, that 
for a polarization which is not in the plane one has to 
compute the actual field direction by use of the Fresnel 
formulae ^3|. The value of the dielectric constant, how- 
ever, which has to be chosen for such a calculation, is 
very close to 1. Tables I and II also give the optimal an- 
gles Qe and $£• for observing (f , f ) and (tt, 0) in higher 
Brillouin zones. Obviously, by choosing the right zone 
and the right polarization a considerable enhancement of 
the intensity can be achieved. 



V. SURFACE RESONANCES AND THE ENERGY 
DEPENDENCE OF THE INTENSITY 

We have seen in the preceding sections that the orbitals 
from which the ZRS is built emit photoelectrons predom- 
inantly at grazing angles with respect to the Cu02 plane 
(see Figure ||) . In most cases electrons emitted with mo- 
mentum (fcy + Gii, k'j_) will be simply 'lost' because they 
will not reach a detector positioned to collect electrons 
with momentum (fey , fcj_) (it might happen that the pho- 
toelectron 'gets rid of its Gm ' by Bragg-scattering at the 
BiO surface layer - a possibility which we are neglecting 
here). It may happen, however, that the motion with 
momentum k\\ + Gy parallel to the surface 'consumes so 
much energy' that there is hardly any energy left for the 
motion perpendicular to the surface. This will happen 
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whenever the kinetic energy T is just above the so-called 
emergence condition for the reciprocal lattice vector Gu: 



T 



2m 



(fen + Gil 



(25) 



In other words: the kinetic energy is 'just about suffi- 
cient' to create a free electron state with 3D momentum 
(fc|| + G, 0). In this case the energy available for motion 
perpendicular to the surface is not sufficient for the pho- 
toelectron to surmount the energy barrier at the surface 
and it will be reflected back in to the solid. 
On the other hand, if the 'perpendicular energy' is below 
the energy of the Hubbard gap as well, the electron can- 
not re-enter the solid either, because there are no single- 
particle states available in the solid with the proper en- 
ergy. The situation thus is quite analogous to the case 
of the so-called Shockley state seen on the Cu surface: a 
combination of surface potential and a gap in the single- 
particle DOS causes the electron to be trapped at the 
surface of the solid. For an extensive review of the prop- 
erties of such so-called surface resonances see Ref. Q . 
The existence of such surface resonances also in 
a cuprate-related material has been established by 
Pothuizen J25]j . In an extensive EELS study of 
Sr2CuCl202, Pothuizen could in fact identify not just 
one, but a total of 3 such states (in some cases with 
clearly resolved dispersion) in the energy window — 
30 eV. His results show very clearly that such sur- 
face states do exist in Sr 2 CuCl202 - whether this holds 
true for other cuprate materials has not been established 
yet, but for the moment we will take their existence for 
granted. 

To proceed in at least a semiquantitative way we use a 
very simple approximation [ p4| . We decompose the po- 
tential felt by an electron at the surface as 



V{r) = V av (z) + V 1 (r), 

V av (z) = — I dx dy V(x,y,z) 
a Jo Jo 

V 1 (r) = V(r)-V av (z), 



(26) 
(27) 
(28) 



(a denotes the planar lattice constant) and assume that 
Vi can be treated as a perturbation. The eigenstates of 
an electron in the potential V av {z) can be factorizcd: 



with corresponding energy 



E(G\\) = E, 



h 2 {k\\ + G\\ 
2m 



A surface resonance would correspond to a ^(z) which 
is localized around the surface. 

Let us now consider the subspace of eigenstates of V av (z) 
which comprises 



a) the (discrete) surface resonance state ^G n ,ti an d 

b) the continuum of states \&(fei|,fc_i_) which evolve into 
plane waves with momentum (fen , k±) asymptotically far 
away from the surface. The periodic surface potential 
(i.e. the term Vi) then provides a mechanism for mixing 
^G||,/i and the continuum. Moreover, dipole transitions 
of an electron from the valence band into both a contin- 
uum state and into the surface resonance (whereby the 
latter alternative is possibly even more probable due to 
the special radiation characteristics of a ZRS) are possi- 
ble. We thus have all ingredients for the standard Fano- 
type resonance (see Figure |l^) and we would thus expect 
a pronounced peak 
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FIG. 12. Energy level diagram for ARPES via a surface 
state. 

in an intcnsity-versus-photon-energy curve whenever 
the kinetic energy of the photoelectrons approximately 
matches the energy of a surface resonance. Neglecting 
the (presumably small) energy this will happen when 
the kinetic energy T obeys the emergence condition ( |25| ) 
for some Gii. 

A strong oscillation of the intensity of the ZRS-derived 
first ionization states at (tt/2,tt/2) and (0.77r,0) in 
Sr2CuCi202 as a function of photon energy has indeed 
been found recently by Diirr at al. fill . Their data show 
a total of four maxima of intensity in the photon-energy 
range 10 — 70 eV. Much unlike the absorption cross sec- 
tion oscillations seen in EXAFS, these maxima are sep- 
arated by near-zeroes of the intensity for certain photon 
energies. 

One interpretation which might come to mind would be 
interference between the directly emitted electron wave 
and partial waves which have been reflected from neigh- 
boring atoms, similar as the variations in absorption cross 
section seen in EXAFS. If we neglect the energy depen- 
dence of the scattering phase shift at the hypothetical 
scatterers, the difference d in pathlength between the in- 
terfering partial waves would then obey 

d(C+T - K ax ) = 2 ^ 



ff 



where k r „ 



is the free electron wavevector at the 



,th 



maximum of intensity. From the four maxima observed 
by Diirr at al. at (0.77T, 0), one obtains three estimates for 
d: 10.6, 12.6, 9.5 A. This is in any case much longer than 
the Cu-0 bond length of 2 A. The only 'natural length' 
in Sr2CuCl202 which would give a comparable distance is 
the distance between the two inequivalent CuC>2 planes, 
which is 7 A. One might therefore be tempted to explain 
the oscillation as being due to electrons from the top- 
most Cu02 plane being reflected at the first Cu02 plane 
below. This would give a difference in pathway of ap- 
proximately 14 A, the discrepancy with the values given 
above could possibly be explained by an energy depen- 
dence of the phase shift upon reflection. On the other 
hand it is quite obvious that the reflected wave, having 
to travel a total of 14 A through the interior of the solid 
and having undergone one reflection, would have a con- 
siderably smaller amplitude than the primary wave. In 
this picture it is then hard to explain why the intensity 
at the minima is so close to zero. We therefore believe 
that interference is the correct explanation of the strong 
oscillations. 

Let us assume on the other hand, that the maxima are 
due to Fano-resonances of the type discussed above. As 
already mentioned above, we would expect a surface res- 
onance state to form whenever the emergence condition 
( p5| ) is fulfilled approximately for some reciprocal lattice 
vector Gii . 



probability for the dipole transition from the Cu02 plane 
into the surface resonance state is proportional the 'tran- 
sition matrix element' uizrs into a plane wave state with 
momentum (fen +Gm , k± = 0). Neglecting the form of the 
Fano lineshape (which would be hard to compute anyway 
because we do not know the Fourier coefficients of the po- 
tential Vi) we then expect an energy dependence of the 
intensity to be roughly given by 



Gii 



I(hu) « £ \m ZR s(k\\ + G||)| 2 ^(T - ^(fc|, + G,|) 2 ). 

(29) 



Actually injection into the surface state should be possi- 
ble if the kinetic energy is in a narrow window above the 
emergence condition, whence for simplicity we replace the 
S- functions by Lorentzians. Figure |l3| then compares this 
simple estimate to the data of Diirr at al. |LlJ . Thereby 
we have again assumed that T = hv — 8 eV. The agree- 
ment with experiment is satisfactory, given the rather 
crude nature of our estimate for the intensity. Given the 
fact that surface resonances in Sr2CuCi202 have been 
established conclusively by the EELS-work of Pothuizcn 
p5[| we believe that the present interpretation of the in- 
tensity variations is the most plausible one. The ques- 
tion to whether such surface resonance states exist also 
in metallic Bi2212 must be clarified by experiment. 
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FIG. 13. Intensities of the first ionization state at (■£ j) 
(top) and (0.7-7T, 0) (bottom) in S^CuCbCh as a function of 
photon energy. The dots are experimental data from Ref. jnj , 
the line is calculated using (pj|). 

As a rough approximation, we moreover assume that the 



VI. APPARENT FERMI SURFACES 



Let us now consider what happens in the above picture 
if we vary the in-plane momentum fey. It is important 
to notice from the outset, that the group velocity of the 
surface state, vy = ^-(fcy + Gy), is much larger than that 
of the valence band. Thus, if we shift the valence band 
'upwards' by the photon energy hv, this replica and the 
surface state my intersect the surface resonance disper- 
sion at some momentum (see Figure 14). Note that here 
the binding energy of the 02p level must be incorporated 
into the valence band energy. For each fey we now assume 
that the Fano-type interference between the surface reso- 
nance state and the continuum occurs. The maximum of 
the resonance curve thereby will roughly follow the dis- 
persion of the surface state. In an ARPES experiment we 
then obviously probe the intensity along the Fano-curve 
at an energy which corresponds to the shifted valence 
band (see Figure |l4|). Obviously this leads to a dras- 
tic variation of the observed photoelectron intensity: at 
the fen labeled 1, the point where the resonance curve is 
probed is on the ascending side of the resonance, whence 
we observe a moderately large 
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FIG. 14. Creation of an 'apparent Fermi surface' by a com- 
bination of Fano resonance and surface state dispersion. 

intensity. At the second momentum, we are right at the 
maximum of the respective Fano-resonance, so the inten- 
sity will be high. At the momentum labeled 3, however, 
we are already at the minimum of the Fano-resonance, 
where the intensity is approximately zero. We thus would 
observe a dramatic drop in the observed intensity over a 
relatively small distance Ak in the fey-plane. More pre- 
cisely, the sharpness of the drop is determined by the 
width W of the Fano-curve and the difference Av g in 
group velocity between the valence band and the surface 
resonance state: Ak w W/Av g . 

To be more quantitative, we present a simple model cal- 
culation of the ARPES spectra to be expected in the pres- 
ence of a free-electron-like surface resonance. We label 
the surface resonance state as |0) and the continuum of 
LEED states by their 'perpendicular' momentum: 
The Hamiltonian then reads 



H Q = |0)£(fc,|,G|| 



fc, V L 



El* 



H.c), 



h 2 k 2 



(30) 



E{k h Gn) = 



ft 2 (fc|| +G\\ 
2m 



(31) 



We assume that the vacuum consists of a volume spanned 
by N x N planar unit cells and thickness L perpendicular 
to the surface. We write the mixing matrix element as 
V/VT so as to explicitly isolate the scaling of the ma- 
trix elements with L. V should be independent of L and 
since we assume it to be independent of k± we may also 
assume it to be real. 
We then obtain the resolvent operator 



iZo )0 (w) = (0|(a>- l 0+-i/r 1 |0) 
1 

~~ hui - So (fey) - Sp) 



TbU) = Hid fen 

2m 11 




Here vq and uiq have the dimension of energy and fre- 
quency, respectively; only one of them can be chosen 
independently, we choose luq = leV/Ti whence i>o is a 
measure for the strength of the mixing between surface 
resonance and continuum. 

A possible final state obtained by dipole transition of an 
electron from the valence band then is 



I*) = T |0> 



(32) 



where we have again explicitly isolated the scaling of the 
transition matrix elements with L. Both matrix elements 
are assumed to be real. The measured intensity then be- 
comes 

A(uj) = -3 R(lu) 



R(u) = ($|(w-i0+-i/)- 1 |$) 
= Tuj-E(k\\,G\i)-E(Q) 



• 2?E(a) 



(33) 



where we have introduced Ti = Ti/V. Taking into ac- 
count the dispersion of the initial state as well as its finite 
lifetime (which we describe by a Lorentzian broadening 
F) we approximate the measured spectrum for some mo- 
mentum fcii as 



(^-e(fc,|)) 2 + r 2 



A(u + hi/) 



(34) 



A simulated ARPES spectrum is then shown in Figure 
|i~5| . For simplicity we have chosen e(fcy) to be the simple 
SDW-like dispersion 



e(fey) = — t e ff (cos k x + cos(fc y ))^ 



(35) 



with t e ff = 0.25 eV. It has been assumed that this band 
is completely filled, that means we would not have any 
true Fermi surface at all. Despite this, the ARPES spec- 
tra, which is simulated by using ( |34] ) shows a sharp drop 
in intensity at the intersection of the surface state disper- 
sion E(k»,Gn) shifted downward by the photon energy. 
Here G\\ = (0,2tt). 
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FIG. 15. Spectral intensity computed from (|34|.) The 
spectra are calculated for equidistant momenta along a line 
through the origin (0, 0) (lowermost spectrum) which forms 
an angle of 10° with the (1,0) direction. The quasiparticle 
dispersion e(feii) and the dispersion of the surface resonance, 
shifted downwards by the photon energy hv, are given by 
lines. 



intensity in a window of 10 meV below the Fermi energy 
- a plot which is by now a standard method to discuss 
Fermi surfaces. In Figure fi6| a the Fano-curve A[ui) in 
( p4| ) is replaced by unity, and we see the expected spec- 
tral weight map in which the Fermi surface can be clearly 
identified. 



(71,71) 




(0,0) (7i,0) 

i (JW 



The photon energy in the above example has been cho- 
sen deliberately such that the surface state dispersion 
cuts through the quasiparticle dispersion near the top of 
the latter, so as to produce the impression of a Fermi 
surface. Usually this would happen only by coincidence 
for very few photon energies. At this point is has to 
be remembered, however, that the dispersion of the va- 
lence band in cuprate superconductors seems to show an 
extended region with practically no dispersion in the re- 
gion around (tt, 0). This band portion moreover is ener- 
getically immediately below the chemical potential. The 
almost complete lack of dispersion in this region then 
makes it necessary to rely mostly on 'intensity drops' 
of the valence band in assigning a Fermi surface, and it 
is quite obvious that whenever the surface state disper- 
sion, shifted downward by the photon energy, happens to 
cut through the flat-band region for the photon energy 
in question, the resulting drop in intensity will be almost 
indistinguishable from a true Fermi surface. To illustrate 
this effect we have calculated the intensity within a win- 
dow of 10 meV below Ep for the 'standard' dispersion 
with an extended van-Hove singularity 26 



f 



(cos(fc x ) + cos(fc y )) 



+ t' cos(k x ) cos(ky) 
t" 

+ — (cos(2k x ) + cos(2fc y )) 



(36) 



which (for t = 0.5 eV : t' = 0.15 eV and t" = -0.05 eV) 
produces the well-known '22 eU-Fermi surface', together 
with the 'flat bands' around (tt,0). The constant c in- 
corporates a constant shift due to the binding energy of 
the O 2p orbitals. Figures p~6| then shows the integrated 



(0,0) 




(0,0) 



(71,0) 



FIG. 16. (a) Gray scale plot of the spectral intensity im- 
mediately below Ef, computed from (B4j) with A(lo) = 1. 

(b) Fermi surface (dashed line) and constant energy contour 
of the surface state (full line). 

(c) Spectral intensity in the presence of the surface state, i.e. 
(B4) with the true A(uj). 



The latter is shown in Figure [16|b. In Figure |16|c, on the 
other hand, we have included the Fano curve A{ui) for the 
surface state with G\\ = (— 2-7T, — 27r) (assuming a free 2D- 
electron dispersion with lattice constant 2. .82 Angstrom) 
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and a kinetic energy of T — 33eU. The correspond- 
ing constant energy-contour I \ (fc + G) 2 = 33eV is 
also shown in Figure 0. One can see quite clearly then, 
that the contour attenuates the true Fermi surface on 
its 'backside' and strongly enhances the intensity at its 
'frontside' thus creating the rather perfect impression of 
a Fermi surface arc which intersects the (1,0) direction 
at approximately (0.87T, 0). On the other hand in our 
model the true Fermi surface is the one seen in Figure 
|l6|a. The combination of Fano resonance and flat band 
portion near (ir, 0) thus produces an 'apparent Fermi sur- 
face' which is entirely artificial. 

To compare with experimental data in more detail, we 
recall that upon neglecting the dispersion of the valence 
band state as compared to that of the surface resonance, 
the momenta in the fen-plane where intersections as the 
one in Figure |IJ occur have to obey the equation 

r _ ft 2 (fc|l+G) 2 
2m 



where T denotes the kinetic energy of the ejected elec- 
tron. In other words: the resonant enhancement of the 
ARPES intensity occurs along a constant energy con- 
tour of the surface resonance dispersion. If we stick to 
our free-electron approximation, these are simply circles 
in the fey-plane centered on — G. Figure [I?] then shows 
these contours for two different photon energies, namely 
22 eV and 34 ev. Thereby we have again assumed that 
T = hv — 8 eV . Also shown is the standard 'Fermi 
surface', determined at 22 eV. It is then obvious that 
the '22 eU-Fermi surface' always is close to some por- 
tion of a surface resonance contour. This would imply 
that for this photon energy most potions of the Fermi 
surface could be enhanced - it might also be taken to 
suggest, however, that some Fermi surface portions near 
(it, 0) are actually artificial. It is hard to give a more 
detailed discussion, because even slight deviations from 
the free electron dispersion for the surface resonance will 
have a major impact on the energy contours. For 34 eV, 
on the other hand, the Fermi surface obviously is 'cut 
off' near (tt, 0) by the contour for G = 2n (-2,0). The 
Fermi surface seen at 34 eV then might be obtained by 
following the true Fermi contour (dashed line) until the 
intersection with the (— 2, 0)-contour, and then following 
the latter (where the ARPES intensity is high due to the 
resonance effect and the fact there are states immediately 
below /i due to the flat band around (tt,0)). In this way 
one would obtain a 'Fermi contour' which is quite similar 
to the one actually observed at 34 eV. For even higher 
photon energy there are too many G's which contribute 
their own circles, so that one cannot give a meaningful 
discussion. 



(7C.7C) 




(0,0) (71,0) 

(TC,7C) 




(0,0) (7C.0) 

FIG. 17. (a) The full lines give the constant energy con- 
tours for a free electron-like surface resonance state for kinetic 
energy 22 eV - 8 eV (a) and 34 eV - 8 eV (b). The dashed 
line is the 'Fermi surface' seen in Bi2212. 

The question which of the different Fermi surface por- 
tions observed in Bi2212 are artificial and which are real 
clearly needs a detailed experimental study. A very sim- 
ple way would be to vary the photon energy in small 
steps of (e.g.) leV, in which case one should observe 
a continuous drift of the 'Fermi surface' for those parts 
which are artificial. Unfortunately the present under- 
standing of surface resonances in cuprate materials is 
much too rudimentary to make any theoretical predic- 
tion. Not even knowing any details of their dispersion, 
we have to content ourselves with the very much oversim- 
plified model calculations outlined above to show what 
might happen. Even these simple model calculations do 
show rather clearly, however, that simply plotting inten- 
sity maps near Ep or identifying 'Fermi surface crossings' 
near (jr, 0) by drops of spectral weight is not really an ad- 
equate means to clarify this issue. 

VII. CONCLUSION 

In summary, we have presented a discussion of ARPES 
intensities from the Cu02 plane. Thereby we could ac- 
tually address only a few (but hopefully the most impor- 
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tant) aspects of the problem. 

To begin with, the geometry of the orbitals, which form 
the first ionization states of the C11O2 plane, leads to 
a pronounced anisotropy of the photoelectron current. 
This makes itself felt in a strong dependence of the inten- 
sity on the polarization of the incoming light. This may 
be exploited to enhance the intensity seen in an ARPES 
experiment by choosing appropriate photon polarization 
or going to a higher Brillouin zone. We have presented a 
simple model to guess 'optimized' values of the polariza- 
tion and the Brillouin zone to obtain maximum intensity 
for a given k-point. 

Moreover there is a strong preference for a Zhang-Rice 
singlet to emit photoelectrons at small angles relative to 
the plane. This in turn may lead to a injection of the 
photoelectrons into states located at the surface of the 
sample, so-called surface resonances. Such surface reso- 
nances are well-established in at least one cuprate-related 
material, namely S^CuC^OV We have proposed to ex- 
plain the pronounced photon-energy dependence of the 
ARPES intensities in this material by these surface res- 
onances. Next, we have presented a simple model cal- 
culation which shows how such surface resonances may 
lead to apparent Fermi surfaces in the flat-band region 
around (tt, 0). This may be one explanation for the ap- 
parently different Fermi surface topology seen in Bi2212 
at 22 eV and at 34 eV, possibly also in an interplay 
with bilayer-splitting |pl 22 1. In any case the existence 
or non-existence of surface-resonances should be clari- 
fied experimentally, since the identification of any major 
Fermi surface portion near (tt, 0) as being 'artificial' amy 
have major implications for our understanding of cuprate 
materials. 
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VIII. APPENDIX I: CALCULATION OF RADIAL 
MATRIX ELEMENTS 

The actual calculation of the radial matrix elements 
Ri^i prove to be one of the most crucial points when mak- 
ing contact with experimental results. Functional forms 
of effective potentials for the CuC"2 compound are not 



1G 



known. A simple approach with a hydrogen-like potential 
seems doubtful, if one attempts cover at least a certain 
amount of physical reality of the photoemission process. 
A good starting point for an at least qualitative deriva- 
tion of the matrix elements are the radial wave func- 
tions that can be obtained from density functional cal- 
culations, which include information about the electron- 
electron exchange correlation and the screening of the 
nuclear charge. The radial matrix elements used in this 
paper arc calculated in the local-density approximation, 
where one solves the effective one electron Schrddinger 
equation 

Q V2 + v eS W) *°>° = M ■ (37) 

The potential t/Jg (r) is a functional of the electron den- 
sity, given by 

v(r) is the full nuclear potential and u([n];r) the direct 
Coulomb-potential fdr'n (r') / |r — r'|. The exchange- 
correlation functional v^ c ([nt, n jj ; r) is a parameterized 
after Ref. p7| . Equation ( |37| ) is solved selfconsistently 
to convergence of the total energy, resulting in a set of 
energy eigenvalues e„ and radial wave functions R n ,i(r), 
which we used to calculate 

Since we assume emitted electrons to be approximated 
by plane waves at the detection point, additional infor- 
mation on the outgoing wave function is needed in form 
of the phase shift compared to the asymptotic behaviour 
of 



defining the interference of wave functions emitted by 
OPx,y and Cud a; 2_ !y 2 orbitals. The phase shift can easily 
be found by comparing the logarithmic derivatives 

[yt R (r) ) /R (r) = (Wr 11 ™™ (r) ) /R * r (r) 
at a sufficiently large radius where 
d 

0~ r v eS « " °" 
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FIG. 18. Radial integrals R;,;/ (top) and Scattering phases 
S U i (bottom) as calculated from density functional theory. 



IX. APPENDIX II 



Here we present in more detail the calculation of the vectors v of interest to us. For the CEF initial states of 
interest to us, such as p x , p y and d x 2_ y 2, the coefficients c am > take the form 

Cam' — 7= \^m' ,u i ^m' v)i 

v2 

(see Table I) whence we find that 

J2 di m Y lm (k°) = E *i» A ^°) 

m fjb= — 1 

A^(k°) = 4= cHl, [i + v- \\ v) ± Yi^k ) c\l, n - v\ I', -v)) (38) 

We now rewrite the scalar product 

1 

E Y *Je)A l fi = e-A l 
n=-i 
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and by use of the property c 1 (Z, — m; I', —m') = c 1 (Z, to; I', to') we obtain the components of A 1 : 



A ltX = ^=[c 1 (l,v-l;l',v) {Yi^tY^^) - c\l,v + ,u) {Yi,v+i T ij,_(„+i)) ], 



A: 



4, 



^= [c 1 ^,!/-!;^,!/) (n,-i±^ > -( I ,-i)) + c 1 (/^ + l;r^) A= (Y^ +1 ±F i; _ (l/+1) ) 



c 



72 



c 



C i^ ; /V) -^(Y^iY,,.,). 



(39) 



Defining R U ,(E) = e iS >R ltl ,(E) we now readily arrive at the following expressions for photocurrent: 



ANhk 



J = 



(40) 



Here the vector v depends on the type of initial orbital. Straightforward algebra then yields the following expressions 
for the vectors v: 



v(p x ) = 



v(Py) 



v{p z ) 



h R21 



-V3d xy {k°) 
-V3d xz (k°) 

-V3d xy (k°) 

£2^+d 3z2 _ r2 (k°) + VSd x 2_ y2 (k ) 

-V3d yz (k°) 

-Vsd xz (k°) \ 
-V3dy Z (k°) 

/ ^ i +Vid 3z2 _Ak°)) 



(41) 



With the exception of the numerical prefactors, most of the above formulas could have been guessed from general 
principles: there are only transitions into s-like and <i-like partial waves (i.e. the dipole selection rule AL = ±1!) and 
acting e.g. with an electric field in x direction onto a p x -orbital produces only s-like, c? x 2_ J) 2-like and rf 3z 2_ r 2-like 
partial waves (which have even parity under reflection by the z — y-plane), whereas acting with a field in y direction 
on p x can only produce the d xy partial wave. We therefore believe that much of the formula remains true even if more 
realistic wave functions for the final states are chosen. 
Similarly we find for the matrix element of the d-likc orbitals: 



v{d 3z 2_ r 2) = i 



( 

-^R 12 p y + yj ^R32hy(z 2 -r 2 ) 



(42) 



v(rf yz ) = 



--J^Rl2Pz — \J^R32.f5z^-3z — -^pfR?,2fz(x2- y 2) 
~^Rl2Py + \f^ R 32f y (5z^-l) 



(43) 



v(4;z) = 



^ -^RviPz - y ^R-32fbz 2 -3z + -JfRz2fz(x 2 -y 2 ) ^ 



V7 

~j E Rl2Px + 



Rd,2f: 



32jxyz 



R32fx(5z 



-1) 



(44) 



18 



v (d x 2_ y 2) = i 



V 



D + 

-^Rl2Py + ^#32.^(5*2-1) + 
-^R32fz(x 2 -y 2 ) 



-3xy 2 



[#32/; 



3yx 2 —y 3 



(45) 



v (d xy ) = i 



14 ^32. fyx 2 



^75^ 12P y ~ ^#32.4(5^-1) + 

■^Rl2Px - -j^R32.fx(bz 2 -l) - \f^R32.fx 3 



V 



75-R32/: 



/ 



(46) 



With the exception of the numerical prefactors, most of the above formulas could have been guessed from general 
principles: there are only transitions into s-like and cZ-likc partial waves (i.e. the dipole selection rule AL = ±1!) and 
acting e.g. with an electric field in x direction onto a p x -orbital produces only s-like, d x 2_ y 2-\ike and d3 Z 2_ r 2-like 
partial waves (which have even parity under reflection by the z — y-plane), whereas acting with a field in y direction 
on p x can only produce the d xy partial wave. We therefore believe that much of the formula remains true even if more 
realistic wave functions for the final states arc chosen. 
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TABLE I. Intensities of the ZRS relative to k = (f, f) in higher Brillouin zones k + G for 22 and 34eV photon energy. 
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TABLE II. 


Intensities of the ZRS relative to k = 


(JT,0) 


m higher Brillouin zones 


k + G for 22 and 34eV photon energy. 
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